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The Effect of the Topography of Surfaces on Adsorption of Gases at Supercritical 
Temperatures on Heterogeneous Solid Surfaces 
Wpływ topografii powierzchni heterogenicznej na adsorpcję gazów 
; w temperaturach superkrytycznych 


BunuaHne Tonorpacbuu SHepreTM1eCKM HEOĄHOPOTHBIX NOBepXHOCTEŃ Ha BEANMIKNY AĄCOPÓLNKU 
B CyliepkpHTWIECKKNX TeMNepaTypax 


INTRODUCTION . 


‘Surface heterogeneity is generally described by an expres- 
sion giving the differential distribution of adsorption sites 
as a function of their adsorption energy. This function, |how=. 
lever, is a satisfactory description of this distribution only 
in baie of an ideal adsorbed phase. In cases where any depar- 
ture from ideal behaviour of the adsorbed.phase exists, the 
differential distribution of adsorption energies is only a 
global description of surface heterogeneity,. $ 

_In such a case another property.of the system which wiil 
influence the behaviour of heterogeneous adsorption is the to- 
pography of solid surfaces. It can be argued that the spatial 
distribution of sites with various adsorption energies over an 
adsorbent surface should influence the global distribution 
function in some way. The role of this topographical distri- 
bution of adsorption sites has not been given much considera- 
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tion in the past, Most publications have assumed a patchwise 
energy distribution [1] with no attempt being made to consider 
the role of the real topography of surfaces studied in physi- 
cal adsorption, 

The alternative topographical model, proposed by Hill [27] 
for the case of localized adsorption, assumes that adsorption 
sites of different energies are distributed randomly over an 
adsorbent surface. This model was subsequently. used by Tompkins 
[3]to analyze the behaviour of a two-dimensional van -der Waals 
gas on heterogeneous surfaces.) : 

Steele [4] was the first to SAY an attempt at discrimina= 
ting between these two topographical models using a virial des- 
cription, but some details of his computations were criticized 
by Adamson [5]. 


Obviously, a real solid surface always exhibits some moder=+ 


ate topographicalicorrelation of adsorption sites and lies be- 
tween the twojextremes described above. Recently, Rudziński [6] 
and Rippa and Zgrablich £72] have studied this spatial correla-. 
tion using virial formalism. This approach is limited to the 
region of low coverage and no investigation has been reported 
of high-coverage adsorption on surfaces with a moderate degree 
correlation of sites of equal adsorption energies, 

\The purpose of this work is to investigate the role of sur- 
face topography for a surface with a moderate spatial correla- 
tion of sites of equal. adsorption ‘energies at the region of 


medium coverage. > 
` 


I. TOPOGRAPHICAL DISTRIBUTION OF SITES 


The overall adsorption isotherm, v (p) for a patchwise mo- 
del of heterogeneous surfaces is given by 


v (p)= Je(p.£) xee | 1.1 
& 


where @(p,&) is the local isotherm on a surface patch having 
adsorption energy equal to €. X(€) is the differential dis- 


tribution of adsorption energy, and dh is the range of possi- 
ble variations of & . 
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Various analytical approximations for X(€) have been propo- 
sed [8-16]. We can argue that the-topography of the surface 
has no effect on the adsorption process until the mutual in= 
teractions between adsorbed molecules start to play a signifi- 
cant role, This is in range of relatively moderate and higher 
coverages of surface. 

In this range we define a new function X (LE) as the \differ= 
jential distribution of the nearest neighbouring sites among 
various adsorption energies T with respect to an adsorption 
site with an adsorption energy © . X%(T,€)\can then be looked at 
as a local distribution. In case of the ideally patchwise to- 
pographical distribution, (U €) is just a Dirac delta fun- 
oeeo 


EN ZUTE) -8(TE) ) 4.2. 
For the fully random topographical distribution, X (TE) should 
be independent of AE Thus 

x(T,€) =$(L) RJ „AD 
Ve are interested here in an intermediate case for which we 
propose the following function 


x (t,€) = 5 ię 14 


"This is a bell-shaped function centered at U= with a ae 
of the order of p - In general, P will -be dependent on E , but 
here we only consider the simple case where p is zada 
of E; 

The potential of an average force acting on an admolecule 
from its nearest neighbouring admolecules is given by ákT Afp). 
whore A(p) denotes the relative coverage of the nearest neigh- 
bours sites of a given admolecule. However, A will, in general, 
differ from both v andQ . In a first-order approximationA €, p) 
can be represented by 
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© s LE 
A(£,p) Ne) | elen) elt) 


E 


where N(@)= 1 + exp dy Jis a normalization coefficient of 
the local distribution. function (1. 4). Using the same expan = 
sion for Fermi-Dirac integrals we can evaluate A E, p) with any ` 
desired accuracy. For thè case of an intermediate topograph~ 
ical distribution, it should be sufficient to retain ey the 
first three terms in this expansion. Hence, 


A(E,p)=8(E,p) + 2(tn2)p aBlpe) + ¥ p* erig 1.6 


where (2) is the local isotherm. In this work we use the Fowler= 
-Guggenhéim isotherm as given below: 


o(p.e)- fa + po emt" 


Equation (1.6) is ‘ebteined dduk that N (€)= 2. This last 
approximation was accepted because we want the truncated se~( 
ries (1.9) to fulfill the necessary condition) 


ahora a oor T 
p-0 zy 


The equation for the overall adsorption isotherm „(M) (p) can 
now be written in the following form: 


vp- f* om KE KK oi p (cote Na 


eat 


Equation (1.9) can be integrated by parts 
vti) (P) = vı (P) + V2(P) i 90443530 


where 
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w (P) = h $ = exp ZE zaj „i JL. 


10) fm (Boo. 


|We shall introduce the Jati transformation | 


and 


ba $ akT AlE) ot SCE, 


"> =t ; 

£ exp Ee exp (4 
res 4kT. (ne -£))¢-0 
"Thus, © takes now the form of a Langmuir equation 


Ə (t p) = (2 + exist sj 


In terms of the new variable t, dint has the form of a 
-Dirac integral ? 


: exp (Soe =) 
er fee ET Ć 


where `ò 


o(t je) = exp |- ett + to - act A) 


Expanding Ọ itang t =JL , we obtain 


ble, B= A a i AEA ee 


1.12 


1.13 


1.14 


1.15 


1.16 


Fermi- 
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E SET the new variable x= a3 we can write. 


Val) = 5, o Piel ya G n (1+ CO 


We shall now apply a further B AE to the integrals 
appearing in Equation (1. 20). 


sarin to Equations (1.14) and (i. 5) we have x 
8 OE) _ 1.21 
=n on ner -0G - . 


Furthermore, we arolinterested in the region of moderate and 
higher coverages of the surface, where interactions between 
admolecules play an important role and the effects of surface 
topography ‘become important. Thus, when evaluating the inte - 
grals in Equation (1.20), we shall accept that in this region ` 
of overall surface coverages, the sites with the smallest ad=- ` 
sorption energy are half-covered. The integration limits are 
then from zero to plus jebnicy; | and we have [18] . 


co s aż 
x dx ; 
= n!B - 1.22 
(1+e*) (i+e"*) z 
where B, is expressed by Riemann 's zeta functions Z(n) 
1- . 
B, = (1-27")z(n) 1.23 


Ce, BL = 1/2, in2, SĘ Hisa; for n = 0,1,2... respectively. 
This approximation imnediately limits the range of applica- 
tion to moderate coverage, the region where the main interest 
of this publication lies. For Q = 0.5, the surface is hali cov= 
ered lat zero energy. The change of integration limit involves 
no error, However, for @ = 0.8 the integral in Equation 1.20 
gives roughly about 0.29, 0.20 and 0.10 between y = 0.0 and 
1.4 for n = 0,1 and 2, respectively. These values, which re- 
present the errors introduced by the change of limit are con- 
pared with the total values (integration limits from 0.to eo) 
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of 0,50, 0.69 and 1.64, Furthermore, because of the symmetrical 
nature of the integrals for n = O and 2, the errors introduced 
by the limit change for O = 0.2 will be of the same magnitude 
as those for Q = 0.8. Thus, a significant error may result if 
is greater than 0.7 or less than 0.3. ~ 
Retaining only the first three terms of Equation 1.23 we 
obtain - : 


wto [1/2 oce) +kTlm2 a oer S Ft) ZY | zy 


t=u 
Combining Equation (4.11) and (1.24) we |finally get 


> (p) [Peo] | +$ cP)" exp|2ér (1 +5 ej] 
* 4-2 (m2) r (1-3) | +e A GE (m +1 -3)]+ 


i SIE 
ir [rA-5)] h +3 ge zk] ! 1.25 


In the limit p-20, y(ltjis the Dirac delta function § (T-£) 
and Equation -(1. 25) reduces correctly to the isotherm equation 
of ideally patchwise surfaces. 


RSE CAI A ERs Z: 
ae -2(iżyri -5) + Er al 


Furthermore,in the Limit $ = Tc/T > O the isotherm of vW or 
reduces to equation 


v= E + | i (e,| > pay? Bor" 1.27 


which is the exact form of v. (p). obtained for Langmuir local 
isotherm. Thus, when the mutual interaction between /adnole- 
cules vanishes, the topography of surface has no effect on ad =- 
sorption ‘tea, 

If there is no spatial correlation between adsorption sites 
of equal adsorption energy, the distribution of adsorption|sites | 
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among different adsorption energies is exactly the same as 
for the whole heterogeneous surface. With the superscript (R) 
“denoting the random Kojak] distribution „of adsorption 
9070008 the local isotherm R (E; ee bg; the form. 


g*) €) - [2 + Ke exp “ese pun hae 1.28 


and. the corresponding overall isotherm v®) (p18, then, given 


s Xp) = [oe] ae SĘ H. 


D-2 tmar +P] 


1.29 


3 z. ~ i y 
"II. THE ISOSTERIC HEAT OF ADSORPTION 


Since the total coverage v is a function of p and T there- 
fore, at constant v, we shall write the equation for the iso- 
steric heat of cag ck Qt in the following form ! 


Pouch: 0/207) 
{av EE 


DE p.T)= k 2.1 


Using equation (2. 1) we may dlcdiace the isosteric heat 

adsorption for the three theoretical isotherns, vP, vt) MO 

developed for the different topographical models of surface. 
(a) Isosteric heat of PROW on ideally patchwise hete- 


rogeneous surfaces Q, R ie 


Qa: (P) (P.T) s kT [e RECO + ah hok in T5)+ 


: + e?e (p, ( p) SE. - E -6-5)]] / [e (1-0) 4 r(v on | 
where H à 


06) „5 TO , ate "2.3 
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n= 1 - 256 WC ZK aw, | 2.4 


e- B® Pe] o | 2.5 


(b) Isosteric heat of adsorption on fully random heteroge- 
neous surfaces Q., RI 16 


o Pp. T) » ee (e (1-0) o®), ię: a) (a Ar ln res) 


+ jo 4r5, sA (Ee y (21n2 - g? J Ve (1-0) + 2.6 


. and 


+r (v® - e) | 2.6 
where. n 4 

ghg, 69, sę vi 27 

e's [e R) e9] 2.8 


(c) Isosteric heat of adsorption oniheterogeneous surfaces 


exhibiting moderate topographical correlation of sites ee is 


a, P.T) = kT l: 1 + r(y (W_ oht- ln 2. + 


3 2 (inz) & ey): $ oxp | 26 r (> z). (ey i 

: 3 eÉ eriy | -Ê (2 +5) - ))- 

„82 płoć x (| deb EJ e 
Pzp)||/G E a ea 
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j 6 (1-9) 9 2.10 


sów ` Ç ET 


os . 3 - 
Binp ~ OOF BUG: acy 


te ete oe Ee. £ (p) ią fide. | si 
+ 2e(1-0) BE s - 30 + 1) -3 | ite 

pP ne sift a 
Beas] | - 2» 
and: Fate ŻA 3 
e = [ol .6)]. oC. 


In the limit T ot Q ,(P =O, Do =Q 02 and 26) 
ipo, we have Gee |= a, ¢ Ą 


III. RESULTS AND DISCUSSION 


First, we shall discuss the behaviour of the overall adsorpt# 
ion isotherms; v'P. v ) and v M „ All these isotherms can be 
written in the following general fora 


vl) (p) = SĘ (p. e + + aË) CJA 3.1 


ad s 
- 


where i = P, R, M, and ao) is, in general, some function of 
(1/1. r and (PART). In other words, every overall isotherm 
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equation is a linear combination of the appropriate local iso- 
therm for a homogeneous surface with the most probable value 
of adsorption eneray* ES and a Freundlich isotherm function 
of the form Bo) . 

. The. last lern h (2 (y in Equation (3.1) illustrates the 
limitations which arise from the several approximations accep- 
ted here. It is, clear that our equations can only be applied 
in the region of moderate surface coverages, in which region 
the Freundlich isotherm can be applied in the absence of inter- 
actions between adsorbed molecules, Thus, although we present 
our numerical results over the whole region of surface coverages, 
we should remember that at very low and very high coverages our 
results are not quantitative, . 

We can now discuss the role of surface topography by compar= 


_ ring the overall adsorption isotherms’ for the two extreme cases 


of surface topography: the patchwise, and the random topography 
of’ adsorption sites. Furthermore, we shall confine ourselves to 
‘the region of supercritical temperatures. 

Taking (r, ft) = 0.5, we have evaluated’ v(P) (p) and v2 (p), 
for three different values of the dimensionless heterogeneity 


'parameter'r: r = 0.01, r = 0.1 and r = 0.3, The results are 


showr in Fig. 1. In Fig. 2 another comparison is made between 
these curves, evaluated ‘now at the constant value r = 0.1 and 


various values of (T/T): (T/T) = 0.1, (T/T) = 0.5 and 


(Te/T) = 0.75. 

A general conclusion, which can be drawn from these two fige 
‘ures, is as follows: surfaces with random topography of ‘sites 
behave more like homogeneous surfaces than surfaces character= 


| ized by patchwise topography. In other words, the patchwise to- 


pography exaggerates the effects of surface heterogeneity. It 
also appears that at low surface coverages, adsorption on sur-. 
faces with random topography is smaller than on those naving a 
patchwise topography. This relationship becomes reversed at 
high overall coverages of surface. The explanation for this ob- 
servation is as follows: 

For high values of, in general, 8(p) > v(P)- Thus, at 
small adsorbate pressures when adsorption is occurring mainly 


in case er the POZOSTA distribution funętion, the 'low= 
est energy Ẹ o is simu taneously the most probable energy, 


The Effect of the Topography of Surfaces... > 116 


Relative surface coverage 


Adsorbate pressure P/K* à 


- Fig, 1. Calculated overall isotherms; 1 - vP) p), 2 - v® oo), 
at (T/T) = 0.5 for three values of the dimensionless hetero- 
geneity parameter r" A 


rage 
x=. 


Relative surface cove 


Adsorbate pressure P/K' 


Fig. 2. Calculated overall isotherms; 4 - v(P)p and 2 ~ v&®Xpy, 
evaluated for r = 0.1, and for three values of (T/T) 


on highly energetic sites, adsorption on surfaces with random 
topography will be smaller than on those with patchwise topo- 
graphy. An opposite effect will be observed at high coverages, 
where adsorption will occur mainly on low-energy sites where 
_e(p) < vip). ] (rt 
- Figure 3 shows some limitations of Eq. (2.25), arising from 
cutting the expansion (1.6) after the third term. Values of 
(Per) higher than 0.4 could not be accepted since MO) (p) 
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Relahve surface coverage 


Qi Q2 
Adsorbate pressure P/K" 
Fig. 3. Calculated overall isotherms v™) p), evaluated at 


r = 0.01 and (T_/T) = 0.5 for three different values of(Ẹ/RT); 
in case of(P/RT = =) the curve was evaluated from eq. 2.2 


` 


Relative surface coverage. 


01 Q2 03 04 


Adsorbate pressure P/K' 


Fig. 4. Calculated overall isotherms vt) (5) at (/RT) = O.1 
and (T/T) = 0.5 for three values ofr” 


would then exceed the values of v(R)(5) at higher adsorbate 
pressures. Obviously, the highest acceptable value of (P/RT) 
depends on both r and (T/T). This limitation can be removed 
as far as necessary by taking higher terms of the expansion 
(1.6) into account. Fig. 4 shows how the overall: surface het= 
erogencity: affects the shape of the overall adsorption iso- 
therm in systems with some moderate spatial correlation of 
sites. The next figure, Fig. 5, provides a similar illustra- 
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Relative surface coverage 


Adsorbate pressure Pjk’ È 


- Fig. 5. Galculated wórsii isotherms v(" ) (p) a at (Pyar = 0. i 
and r = 0.01 for three different “values of ( 


+ 060-060 0 
Relative surface coverage VP) 
Fig. 6. The comparison of theoretical, isosteric heats of ad- 
sorption between surfaces with patchwiee and Sandoli dz ść ; 
(8) 


of sites; ail curves were evaluated at T = 100°K, - ua 0. 
and C/T) = 0.5 for three values of. opa 


tion of the effect when the value of (1/7) dings; 

Figures 6-10 show the effect of surface topography in the 
isostericheats of adsorption. All these curves have been eval-- 
uated for| the temperature T = 100°K, The reason for these 
figures is that isosteric heats of adsorption are probably the 
most extensively investigated adsorption characteristics next 


to the isotherms themselves. In various publications, heats of 
adsorption are commonly presented as a function of surface co- 
verage in the way we mz: them here, 


"2 
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Figure 6 shows the comparison of the isosteric heats of ad- 
sorption for the two extreme cases of surface topograpny: for 
the patcnwise, and the random topography of sites. Curves such 
és shown here are often observed experimentally [13-23], and 
their interpretation is as foiiows: the first, dyucreasing 
branch of these curves, corresponds to the coverage region in 
which the effect of decreasing adsorption energies prevails 
over the effect of increasing energies of interaction between 


kcal 
"mole 
5 


St 


Heat of adsorption Q 


Relative surface coverage VP) 


Fig. 7. The comparison of theoretical isosteric neatś of ad- 
sorption between surfaces with patghwise and random topography; 
all curves were evaluated at T =100 K, bo =.0.0, and r=Q,1 

for three different values of (T/T) 


; kcal 
Heat of adsorption Qst mole 
3 


8 


04 06 08 10 


Relative surface coverage V(P) 
Fig. 8. Calculated isosteric heats of adsorption for surfaces 
of moderate spatigl correlation of sites; all curves were e 
‘uated at T = 100 K, = 0.0, r = 0.01, and (T/T) = 0.5 for 
three 8ifferent values of (PRT) 
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Heat of adsorption Qa SH 


04 06 08 10 
-Relotive surface coverage VP) 


Fig. 9, Calculated isosteric heats of adsorption for surfaces 

with moderate page! correlation of sites. All curves 

evaluated at T = 200K „ E = 0.0, (T_/T) = 0.5, and Crapo. as 
for three different Selues of r 


$ 04 06 
ię Relative surface coverage V(P) 


Fig. 10. Calculated isosteric heats of adsorption for surfaces 
' with moderate a ti correlation of sites, All curves. were 
evaluated at T = 100°K = 0.0, r= 0.01, and ng = Oot 
for three different values of /1) 


“adsorbed molecules. The right-hand side, ORRETAN branch of 


these curves, bocomes linear as predicted by the Fowler=Guggen= 
heim equation for the whole region at supercritical a h 


res and on homogeneous surfaces. 

(Figure 7 shows that in some cases only the decreasing branch 
of the isosteric heat of adsorption will be found. Sinco such 
curves are sometimes used as an approximate estimate of the 
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overall distribution of adsorption energy we see how risky 
this procedure can be. ] 

Figures 8-10 illustrate the behaviour of the isosteric heats 
of adsorption in various physical situations, in systems with 
some mediate topographical correlation of sites. 


CONCLUSION 


, 


By introducing some approximations, analytical expressions 
have been obtained for both the overall adsorption isotherms 
and for the overall isosteric heats of adsorption; in adsorp- 
tion at supercritical temperatures for three basic topographi- 
cal models of heterogeneous surfaces: surfaces with patchwise 
tepography of sites, surfaces with_spatial correlation of sites 
of equal adsorption. energies, and for surfaces with a ran- 
„dom topography of sites, zy 

He find that in both the overall adsorption isotherms and 
the theoretical heats of adsorption the effects of surface jhet- 
jerogeneity are most pronounced in case of surfaces with a 
patchwise topography of sites, : 
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STRESZCZENIE 


Rozważono konsekwencje przyjęcia trzech różnych modeli topo- 
peatas powierzchni (płatowego, random i pośredniego) w zjawis= 
ach adsorpcji na powierzchniach heterogenicznych w obszarze 
temperatur asg LE (pi ay Badania przeprowadzono w oparciu 
o tzw. całkowy formalizm opisu izotermy adsorpcji globalnej, 
przyjmując do ba geta adsorpcji lokalnej model Fowlera-Guggenhei- 
ma. Przy pewnych uproszczeniach o charakterze matematycznym o= 
trzymane zostały analityczne wzory dla izoterm adsorpcji i cał= 
kowitych izosterycznych ciepeł adsorpcji. Stwierdzono, że na po- 
wyższe wielkości termodynamiczne najwiekszy wpływ ma niejedno= - 
rodność powierzchni o topografii płatowej. : 
` $ ij z» 
, 48 Peswue 7 
_ B zanno paóore paccyx46HH pe3YIBTATH upuMeHeiMA TpeX pa3HLX 
monene Tonorpaduu noBepXHOCTU (aockyTHo%, cnygażHoii u cepeZuH— 
Holl) B TEOPETMYCCKUX WCCJIEĄOBAHUAX A1COPÓNUM HA BSHEpreTuyecKu— 
-HGOAHOPOAHNX NOBEPXHOCTAX. HCCJIENOBAHKA MpOBeAeCHH Ha OCHOBA- 
` HUM MHTEPPSAbHOTO QOPMAJIM3KA OMUCAHKA M3OTEPMH oome ti ancopouuu. 
B KayveoTBe JOKAJIBHOŃ W30TepMH AĄCOPÓNUM TpuHATO ypaBHeHue QOB- 
nepa-TyrreHxeńwa. [Ipu H6KOTOPHX MATEMATKUECKAX JNPONEHUAX NONY- 
U8HO AHANKTUUECKKE QODPMYJK JIA UsOTEpM ĄCOPÓNKM U OÓNUX H30- 
cTepuyeckux TENIOT AĄCOPÓNKK. JCTAHOBJIEHO, YTO camoe ÓOJBMOG 
BIWAHMC HA STU TEpMOĄKHAMKUECKKE BEJKUWHHN UMEET, HEOAHOPOTHOCTS 


NOBepXHOCTM C JIOCKyTHOŃ ronorpadneli. 


